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ABSTRACT 

We use a deep XMM-Newton Reflection Grating Spectrometer observation to examine the X- 
ray emission from the core of the Centaurus cluster of galaxies. We clearly detect Fe XVII 
emission at four separate wavelengths, indicating the presence of cool X-ray emitting gas in 
the core of the cluster. Fe ions from Fe XVII to XXIV are observed. The ratio of the Fe XVII 
17. 1 A lines to 15. OA line and limits on O VII emission indicate a lowest detected temperature 
in the emitting region of 0.3 to 0.45 keV (3.5 to 5.2 x lO*" K). The cluster also exhibits strong 
N VII emission, making it apparent that the N abundance is supersolar in its very central 
regions. Comparison of the strength of the Fe XVII lines with a Solar metallicity cooling 
flow model in the inner 17 kpc radius gives mass deposition rates in the absence of heating 
of 1.6 — 3 M0 yr^' . Spectral fitting implies an upper limit of 0.8 M0 yr^^ below 0.4 keV, 
4Mq yr^^ below 0.8 keV and 8 M0 yr^^ below 1.6 keV. The cluster contains X-ray emitting 
gas over at least the range of 0.35 to 3.7 keV, a factor of more than 10 in temperature. We 
find that the best fitting metallicity of the cooler components is smaller than the hotter ones, 
confirming that the apparent metallicity does decline within the inner 1 arcmin radius. 

Key words: X-rays: galaxies — galaxies: clusters: individual: Centaurus — intergalactic 
medium — cooling flows 



1 INTRODUCTION 

In the centres of many clusters of galaxies the mean radiative cool- 
ing time of the intracluster medium (ICM) is short, typically less 
than I Gyr (e.g. |Voigt &. Fabian|2604) . This is also where the ICM 
is significantly cooler than the outskirts. It is therefore expected 
that a cooling flow l jFabian|1994^ should be formed, where mate- 
rial cools out of the ICM. 

High spectral resolution X-ray studies of nearby clusters of 
galaxies using the Reflection Grating Spectrometer (RGS) instru- 
ments on XMM-Newton has revealed a lack of cool X-ray emitting 
gas in these objects ( Tamura et al.|200Ia|b [Peterson et al.|2001| 
IKaastra et a"Ll|200T| |S^elUou et al.||2002| [peterson et al.||200'3l ^ 
There are emission lines observed from gas down to a half or a third 
of the outer temperature, but with very little gas at cooler tempera- 
tures (see Pet erson & Fabian|2006| for a review). The main spectral 
indicator of the lack of cool gas is the weakness of Fe XVII lines. 

A cluster where there has been tentative evidence for the 
existence of cool gas is Abell 2597 ( [Morris & Fabian||2005[ l, 
which shows possible Fe XVII emission in its RGS spectrum. 
NGC 4636 emits strong Fe XVII lines, but these are consistent with 
a small range of temperature as the group is cool ( [Xu et al.|2002| l. 



E-mail: jss@ast.cam.ac.uk 



NGC 5044 is another cool system showing strong Fe XVII and hint- 
ing at O VII emission ( [Tamura et al.|2003) . 

The Centaurus cluster, Abell 3526, is a nearby bright cluster 
of galaxies. It lies at a redshift of 0.0104 ( [Lucey e t al. 1986!), and a 



2-10 keV X-ray luminosity of 2.9 x lO"*^ erg s~"'"( [Edge etal.|1990 1, 



It lies at a relatively low Galactic latitude of 21 .6° with a Hydrogen 
column density towards the cluster of ~ 8.6 x 10^" cm^^ (derived 
from the H I maps of Kalberla et al.[2005| l. 

A cool X-ray emitting gas component has been seen in this 
object by low spectral resolution CCD quality observations. This 
component was first seen by ASCA ( [Fukazawa et al.[[l994[ l and 
ROSAT ( [Allen & Fab ian 1994), who found cool gas with a tem- 
perature of ~ 1 keV in addition to a hot ~ 4 keV component. 
[Ikebe et al.[ jT999j put forward a two-phase model fitting A5CA 
and ROSAT data, claiming 1.4 and 4.4 keV components over the 
inner 3 arcmin. A colour analysis of ROSAT data jSanders et al.[ 
[2000[l gave a me an temperature of 1.7 keV in the central 90 arc- 



AUen et ar] ( [2001| l found an improved fit to their ASCA spec- 



tra with two or more temperatures, with temperatures at 3.2 and 
1 .3 keV. There was also evidence for a hot or powerlaw component 
(seen previously by [Allen et aL|2000[and[Di Matteo et al.|2000l l. 
BeppoSAX observations ( Molendi et a l.|2002^ found a peak temper- 
ature of 4 keV, dropping to below 2 keV in the core, and evidence 
for a hard component in the PDS instrument. 
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The excellent spatial resolution of Chandra provided a more 
complete picture of the temperature of the gas in the central re- 
gions ( [Sanders & Fabian"p 002 l. This observation showed a plume- 
like feature in the core of the cluster, emitting soft X-rays. Two- 
temperature component fits to projected spectra extracted from the 
plume were preferred to single-temperature fits, with temperatures 
of 0.7 and 1.5 keV. Outside this region the temperature increases 
quickly to the east to around 3.7 keV. To the west, there is a plateau 
of cooler gas at around 2.5 keV, before the temperature rises at a 
radius of around 190 arcsec. Deeper Chandra observations jFabian| 
|et al.|2005] l confirmed this picture, highlighting the clear east-west 
asymmetry of the temperature distribution. XMM-Newton EPIC 
data confirmed the ~ 0.7 keV component in the core ( [Sanders "&] 
[Fabia n 2006a l, also finding the temperature of the ICM to drop to 
around 3.4 keV beyond 120 kpc radius. The Chandra observations 
also clearly showed that the metaUicity of the ICM is inversely cor- 
related with its temperature (^S anders & Fabian|2002[ [Fabian et aE] 
[2005|[Sanderr & Fabian 2006a). 

An interesting feature of this cluster is that the metallicity of 
the gas towards the centre is significantly supersolar. The Fe metal- 
licity peaks between 1.5 and 2 Zq ( [Fabian et al.][2005[ [Sanders &[ 
[Fabian|200 6al. Si and S peak around 2Z0, and Ni peaks around 
4 Zq. In addition the metallicity of the gas appears to decline in the 
very central regions ( [Sanders & Fabian|2002[ l. 

The multiple detections of cool X-ray emitting gas in the core 
of Centaurus provides an excellent opportunity to test the picture 
that there is only a range in 2-3 in X-ray gas temperature in clus- 
ters of galaxies. We therefore undertook a deep XMM-Newton RGS 
observation of Centaurus, which, when combined with the existing 
observation, gives a total exposure of around 160 ks. 

In this paper we assume Hq = 70 km s^ ' Mpc^ ' , which gives 
an angular scale of 213 pc per arcsec for Centaurus. We assume 
the Solar metallicities of [Anders & Grevesse[(r989^ . Error bars are 
quoted as 1(7 and limits as 2(J. 



2 DATA PREPARATION 

We processed the two source datasets listed in Table[T]using XMM 
Science Analysis System (SAS) version 7.1.0, with the RGSPROC 
pipeline. We processed the 1st and 2nd order spectra, including 99 
per cent of the point spread function (PSF) and 97 per cent of the 
pulse-height distribution, in order to get as many source photons as 
possible without increasing background too much. We examined 
the observation light-curves for CCD number 9 at absolute values 
of the cross dispersion greater than 1.5 x 10^** (using values of 
FLAG of 8 and 16). The lightcurves were relatively consistent at 
values of around 0. 1 s^ ' , except for some short flaring in the longer 
observation. We filtered this observation, excluding time periods 
with count rates greater than 0.2 s^' . 

Centaurus fills the entire field of view of the RGS instruments. 
We therefore required blank-sky background spectra to subtract in- 
strumental and external backgrounds. As we used a 97 per cent 
pulse-height distribution cut we could not use that standard RGS- 
BKGMODEL tool. Instead we selected five deep RGS observations 
of point-like sources from relatively low Galactic latitude to gener- 
ate background spectra (Table[TJ. We processed these observations 
with the same parameters as Centaurus, cleaning flares in the same 
way, and excluding the inner 90 per cent PSF (where the sources 
lie) to generate the background spectra. 

We combined with RGSCOMBINE the separate background ob- 
servations to make RGS 1 and RGS 2 spectra for the two spectral 



orders. We used RGSCOMBINE to add the spectra and responses 
from the two foreground Centaurus observations, after reprocess- 
ing including the correct attitude values. The background spectra 
contained some invalid spectral channels which did not correspond 
with the foreground spectra, so we marked these channels as in- 
valid in the foreground spectra. We grouped the foreground spectra 
to have a minimum of 25 counts per spectral bin. 

We checked the background spectra were the same by apply- 
ing them to the foreground individually. The backgrounds were in- 
distinguishable in the 6.5 to 27 A range. 



3 SPECTRUM 

We show in Fig. [T| the full fluxed spectrum using the 99 per cent 
PSF with 97 per cent of the pulse-height distribution. A 99 per cent 
PSF corresponds to roughly the inner 160 arcsec width in the cross- 
dispersion direction. This was created by combining the first and 
second order spectra for both observations together with the task 
RGSFLUXER, subtracting the combined background spectra. We do 
not show the spectrum at wavelengths longer than 26A as the back- 
ground becomes increasingly important. We note that the output 
from RGSFLtJXER is for display purposes only. We do not use it to 
obtain quantitative information. 

In the plot we also show the spectrum extracted from the in- 
ner 90 per cent of the PSF (which corresponds to around 60 arc- 
sec width in the cross-dispersion direction). In addition we plot 
the difference between the two spectra, which is equivalent to 
the spectrum extracted between 60 and 160 arcsec in the cross- 
dispersion direction. In the lower panel we show smoothed model 
APEC ( [Smith et al.|200l) l spectra for plasmas with Solar metallic- 
ity at temperatures of 0.5, 0.7 and 1 keV for comparison. These 
models are plotted with arbitrary normalisation to have roughly the 
same range in values. Note that the spectrum also contains lines 
from gas hotter than 1 keV. 

The spectrum shows a variety of emission lines, from N, O, 
Ne, Mg, Si and Fe in several ionisation states. Most interesting are 
those lines indicating cool X-ray emitting gas, particularly strong 
in the 90 per cent PSF spectrum. Fe L lines from Fe XXIV down 
to Fe XVII are observed. The strong N VII line, in particular, is 
interesting. 

We show a zoom-up of the spectrum between 12 to 26 A from 
the 90 per cent PSF in Fig. [2] The plot shows we clearly observe 
several distinct Fe XVII lines. These lines are relatively narrow, in- 
dicating that they come from a relatively small region (See Section 
[3.1[ >. Absent from the spectra is any evidence for O VII emission 
(which would appear at 21.8 and 22. 2A at this redshift), which is a 
strong indicator of gas less than 0.2 keV. 

Different lines are sensitive to gas at different temperatures. 
The strength of lines from a particular ion are governed by how 
much of an element is in the form of that ion, the abundance of the 
element and the density of the gas. We plot in Fig.|3]the fraction of 
an element which is in the form of a particular ion as a function of 
temperature, for those ions from which we observe lines (plus O VII 
which we do not observe). These results are from the calculations 
of ionisation equilibrium of [Mazzotta et al.[ ( [1998) (which are used 
by the APEC model). 

We show in Fig. |4] the image formed by plotting the disper- 
sion angle of detected photons against the cross-dispersion angle. 
Photons outside of the first order spectrum were removed by mak- 
ing an energy-cut using the pulse-invariant CCD energy values. It 
shows the combined RGSl and RGS2 cross-dispersion image for 
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Wavelength (A) 

Figure 1. Top panel: fluxed and combined 1st and 2nd order RGSl and RGS2 data from the two observations. Shown are the spectrum using 95 and 99 
per cent of the RGS PSF. Also plotted is the difference between the two spectra, showing the emission from outside of the very central region. The spectra 
have been binned by a factor of 4 to have 0.04A spectral bins. Bottom panel: smoothed APEC Solar model spectra for gas at 0.5, 0.7 and 1 keV with arbitrary 
normalisation. 



Table 1. XMM-Newton RGS observations of the target (Centaurus) and the five background fields. RA and Dec ai'e the boresight pointing coordinates, length 
specifies the observation length, and the exposure is the RGSl exposure after cleaning for bad time periods (we total 158 and 157 ks exposure on the target in 
the RGSl and 2, respectively). A'h shows the weighted Galactic Hydrogen column density towards the object jKalberla et al.|2005) . 



XMM observation ID 


Target 


RA (J2000) 


Dec (J2000) 


Date 


Length (ks) 


Exposure (ks) 


A'h (lO^o cm 


0046340101 


Centaurus cluster 


12''48""47.93' 


-41°18'43.5" 


2002-01-03 


47.8 


46.5 


8.56 


0406200101 


Centaurus cluster 


12'"48'"43.05» 


-41°18'42.5" 


2006-07-25 


124.3 


111.7 


8.56 


0112210201 


NGC 3783 


ll'"39"'01.70^ 


-37°44'08.8" 


2001-12-17 


137.8 


72.5 


9.91 


0112210501 


NGC 3783 


Ijhggmoi 73s 


-37°44'08.7" 


2001-12-19 


137.8 


124.2 


9.91 


0300430101 


NGC 3256 


10'"27"51.45» 


-43°54'14.0" 


2005-12-06 


134.0 


116.6 


9.14 


0302850201 


MCG-5-23-16 


09''47"H5.43» 


-30°56'57.6" 


2005-12-10 


131.2 


111.7 


8.69 


0305310101 


EC 13471-1258 


13'"49"'52.17^ 


-13°13'36.0" 


2006-01-27 


119.4 


94.6 


5.34 



the two observations, exposure-corrected for bad pixels. The hori- 
zontal axis shows increasing wavelength, while the off-axis angle 
in the cross-dispersion direction is shown vertically. It can be seen 
that some lines, e.g. Fe XVII, are much more centrally concentrated 
than others, e.g. O VIII. 



3.1 Line cross-dispersion profiles 

We can examine this more quantitatively by plotting the profiles 
of individual lines as a function of cross-dispersion angle. As one 
of the most interesting cases, we show in Fig.|5]the profiles of the 
two strongest Fe XVII lines. The 15.0lA line can be strongly reso- 
nantly scattered as its oscillator strength is 2.7. Looking at the top 
panel of Fig.js] which shows the raw profile in a ±0.lA strip either 
side of the line centre, it appears there is evidence for scattering 
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Wavelength (A) 

Figure 2. Top panel: fluxed and combined 1st and 2nd order RGSl and RGS2 data, using 90 per cent of the RGS PSF. The spectrum has been binned by a 
factor of 2 below 20A ( giving O.O2A bins) and 3 above that. Bottom panel: smoothed solar-abundance model APEC spectra at 0.25, 0.5, 0.7 and 1 keV with 
arbitrary normalisation. 




Temperature (keV) 

Figure 3. lonisation equilibrium for lines from observed (solid) and unobserved (dotted) ions as function of temperature in keV ^Mazzotta et al.|1998^ . The 
data are shown here smoothed with cubic splines. 



(a broader ISA profile). However, subtracting nearby continuum 
emission creates mucii more siiarply peaked distributions of emis- 
sion (shown in the bottom panel of Fig.|5j, consistent between the 
two Unes. This is much more in line with the predictions from the 
Chandra Centaurus data of a resonance scattering model ( [Sanders| 
|& Fabian|2006b| l, which predicts at most a few per cent effect (res- 



onant scattering is strong in the very central, cool, dense regions, 
but this cannot scatter radiation to appear to come from a larger 
radius). 

Fitting a Gaussian to the continuum-subtracted profiles gives 
half energy width (HEW) values of 14.2 and 12.3 arcsec for the 
17.06A and IS.OlA lines, respectively. These values are consis- 
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Figure 4. Image of the first order exposure-corrected dispersed spectrum, combining both RGS instruments and observations. The dispersion direction lies 
along the horizontal axis, with increasing wavelength rightwards. The vertical cross dispersion direction shows a one-dimensional image of the central 
~ 5 arcmin of the cluster 
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Figure 5. Line profiles as a function of cross-dispersion angle for the two 
strongest Fe XVII lines. The top and bottom panel shows before and after 
subtracting neighbouring continuum, respectively. 
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tent with the approximate HEW of the RGS mirrors (13.8 and 13.0 
for the RGSl and RGS2, respectively). The 17.06A blended line 
strength appears substantially stronger than the 15.0lA signal. We 
will reexamine this in Section l431 

We show in Fig. |6] cross-dispersion profiles of several other 
strong lines. In each of these cases we plot the profile 0.1 A ei- 
ther side of the centre (except for the broad O VIII and Mg XII 
lines, where we use 0.2A). The general picture is that lines that 
are emitted by cooler gas (e.g. Fe XVII, XVIII, XIX) come from 
more compact regions. One line which differs from this picture is 
the N VII line, which is fairly broad in the cross-dispersion pro- 



Figure 6. Profiles of several prominent lines in the cross-dispersion direc- 
tion. The units are in 10^' counts s"' A^' . Continuum has been subtracted 
from a neighbouring part of the spectrum. 



file. Fig.[T|however shows that most of the contribution to this line 
comes from the inner 90 per cent of the PSF. 
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Table 2. Spectral fit results to the data extracted from the 99 per cent PSF aperture. If two or more results are given, this means the parameter was allowed 
to vary independently in each component, otherwise they were tied together. Upper limits are 2a. For the multi-temperature models the number in brackets 
signify the temperature component number the best fitting parameters apply to. The values of the source sizes are HEW, exclusive of the 13.2 arcsec component 
from the mirrors. 



Model 


VAPEC+VAPEC 


VAPEC+VMCFLOW 
k7;™„ = 0.08 keV 


VAPEC+VMCFLOW 
kT^^jfi free 


5 X VAPEC 


5XVMCFL0W 


kr (keV) 


1.85 ±0.02 


1.86±0.02 




2.04 ±0.02 


3.2, 2.4, 1.6, 0.8, 


3.2-2.4, 2.4-1.6, 




0.766 ±0.008 






0.54±0.01 


0.4 


1.6-0.8, 0.8-0.4, 


jVfi ^lu cm } 












0.4-0.0808 


in 1 -1- n 9 

iU. i It u.z 


10.1 ±0.2 




7-1-09 
y . / It U.Z 


y.KJ It U.Z 


O.O It U.Z 


sizes (arcmin) 


1 1 9 -1- n 0'^ 

i , IZ It U.U-J 


1.19±0.03 




i .^L It U.Ut 


1 Qf\ -\- 00 (\ 9^ 
l.yoitu.uy yV-L)^ 


_ , a / 1 \ n QQ _i_ 




r» Ti \ r\ r\'~i 

0.32 ±0.03 


0.37±0.03 




U.J5 it U.U2 


U.JO It U.UJ {o). 


0.U4 (z), O.iz it 












0.39 ±0.04 (4-5) 


0.06 (3), 0.49 ± 














0.09 (4-5) 


N 


< 0. 18 


< 0.18 




< 0.21 


< U.15 (1-z) 


< 0.18 (1-2) 




^•*-0.8 


0.86±0.12 




2.0it0.3 


1.6±0.2 (3-5) 


3.5 ±0.6 (3-5) 


O 


A n _1_ A AO 

0.32 ±0.02 


0.34 ±0.02 




A /I A _1_ A Al 

U.4U ± U.U3 


A /I -7 _1_ A AO 

U.4/ ± U.02 


A CO _1_ A AO 




0.41 ±0.09 


0.15±0.03 




0.34±0.05 






Ne 


0.45 ±0.03 


0.37±0.02 




0.56 ±0.04 


0.74 ±0.04 


0.80 ±0.05 


Mg 


0.60 ±0.03 


49+0.03 




0.69 ±0.05 


0.91 ±0.05 


1.02 ±0.07 


Si 


1.22 ±0.06 


1.02 ±0.06 




1.34±0.05 


1.60 ±0.09 


1.76±0.12 


Ca 


2.0 ±0.3 


2.2±0.3 




2.1±0.2 


2.6±0.4 


2.8 ±0.5 


Fe 


0.66 ±0.02 


0.67 ±0.02 




0.82±0.03 


0.97 ± 0.02 


1.08 ±0.02 




0.53 ±0.04 


0.29 ±0.03 




0.59 ±0.07 






Ni 


1.81 ±0.08 


1.67 ±0.07 




1.95 ±0.09 


2.34 ±0.11 


2.54±0.15 


norm 


(4.14±0.04) X 10-2 


(3.89±0.07) 


X 10-2 


{3.46±0.06) X 10-2 


See Fig.fTo] 






(2.0±0.1) X 10"^ 












M(M0yr-l) 




7.4±0.5 




8.7±0.5 




See Fig.fni 


Xv 


6435/5068= 1.27 


6494/5069 = 


1.28 


6381/5068 = 1.26 


6362/5069 = 1.26 


6397/5069= 1.26 



4 SPECTRAL FITTING 

The observed spectrum is affected by several factors, listed below. 

(i) The temperature of the observed region affects the ratio of 
the strengths of the emission lines and the continuum. Emission 
lines are increasingly important below --^ 1 keV. Some lines, e.g. 
Fe XVII, are only emitted over certain restricted temperature ranges 
(see Fig. |3|. The spectrum will be the sum of several temperature 
components which means care must be taken in using line ratios as 
temperature diagnostics. 

(ii) The metallicity of a temperature component affects the 
strength of its emission lines. In general the metallicity of each 
element in the plasma can vary as a function of temperature and 
position. 

(iii) The spatial size of the emitting region causes spectral 
broadening. If a temperature component comes from a smaller spa- 
tial region, its lines will be broadened less than a component from 
a large region. Lines are broadened by ^Brinkman et al.|1998| l 



0.124 



(1) 



where A0 is the source extent (HEW; half energy width) in arcmin, 
and m is the spectral order. The XMM RGS response already takes 
account of the broadening for point sources, so this width is on top 
of the ~ 13.2 arcsec HEW PSF of the mirror modules of the RGS. 

(iv) The sensitivity of the RGS is a function of the position on 
the sky. It is less sensitive to components observed from off axis. 

Therefore a complex source like the core of the Centaurus cluster 
would require a significant amount of modelling to account for each 
of these effects. There are also degeneracies when fitting the spec- 
tra. For example, the emission measure of low temperature gas is 
inversely proportional to its metallicity as the continuum is difficult 
to measure. 



Our general technique in this paper is to fit multiple ther- 
mal models to the spectrum extracted from the RGS. Each thermal 
model is smoothed by a variable-width Gaussian (of width in wave- 
length given in equation[TJ to take account of the spatial distribution 
(which assumes the emission region is Gaussian in shape, centred 
on the cluster centre). We tie the abundances of each of the thermal 
components together, unless the data are good enough to constrain 
them in individual components. The thermal components are ab- 
sorbed by Galactic absorption before smoothing with the Gaussian 
model. 

Our modelling assumes that each component has a fixed tem- 
perature, metallicity and spatial scale. In most of the spectral fits we 
make assumptions in our choice of which parameters to vary. For 
example, we do not allow many of the metallicities to vary with 
temperature. 

We fitted the first-order spectra between 6.5 and 27 A, where 
the background is low relative to the foreground. The second-o rder 
spectra were fitted between 7 and 15.5A. XSPEC version 12.3.1 1 Ar- 
|naud| 1996) 1 was used to do the spectral fitting. We use the spectra 
extracted from the 99 per cent PSF to get as complete a picture as 
possible. 

4.1 Single temperature model 

To demonstrate that an isothermal model is inadequate, we show in 
Fig.|7]the best fitting single temperature VAPEC ( Smith et al.|2001| l 
model. In this model the temperature, Galactic column density and 
normalisation are free. The relative normalisations of each RGS 
dataset are allowed to vary relative to the RGSl first order spec- 
trum. The N, O, Ne, Mg, Si, Ca, Fe and Ni metallicities are free 
to vary. The spectrum is smoothed by a variable-width Gaussian to 
account for spatial broadening of the spectral lines (Equation [TJ. 
The 1.6 keV model fails to account for the Fe XVII and N VII lines 
(;K^ = 7474/5075 = 1.47). 



© 0000 RAS, MNRAS 000, 000^000 



Cool gas in the Centaurus cluster 



3 X 10" 



2 X 10" 



10" 




1.6 
1.4 

1.2 

1 




S 0.8 
^ 0.6 




7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 

Wavelength (A) 

Figure 7. Best fitting single-temperature model and ratio of data to model for the 99 per cent PSF spectra. The points are the first and second order data from 
the two RGS detectors and the lines are the best fitting model after folding through the response. The data have been rebinned for display in XSPEC to have a 
signal to noise ratio of 13. The model fails to account for the Fe XVII and N VII lines. 



4.2 Simple two-temperature model 

The most simple extension of a single thermal component model is 
one with two temperature components. This model consists of two 
VAPEC components, each with a variable temperature, normalisa- 
tion, and N, O and Fe metallicities. They are absorbed by the same 
Galactic column density, but are allowed to vary independently in 
spatial size (smoothing according to Equation[TJ. The Ne, Mg, Si, 
Ca and Ni metallicites are tied between the two components and 
allowed to vary independently. We allow the N metallicity to vary 
between the two components because the N VII line is narrow. If the 
hot component is forced to have the same N metallicity as the cool 
component, then the model cannot fit the narrow line because too 
much emission comes from larger spatial regions. Table |2] shows 
the best fitting values of each of the parameters and their uncertain- 
ties. The improvement in is around 1000 over the single tem- 
perature model. The best fitting temperatures, 0.77 and 1.85 keV, 
are very close to the range of values found in the inner arcmin by 
Chandra (Fabian et al. 2005 ). 

The best fitting iron metallicities, however, are substantially 
lower than the 1.5 — 2 Zq peak values from Chandra and XMM 
CCD spectra ( .Sanders & Fabian|2006a| l. In fact all of the metallici- 
ties (except Ca) of each of the elements we measure from the RGS 
spectra are significantly lower than the CCD results. CCD mea- 
surements indicate that the metallicities drops dramatically in the 
innermost region ( Sanders & Fabian 20021, which may coiTespond 
with this low metallicity. We will discuss this further in Section[54] 



The lines from cooler temperature gas are narrower than the 
lines from hotter gas, leading to a larger spatial scale for the hot 
component (1.1 vs 0.3 arcmin). The emission measure of the cooler 
component is only around 5 per cent of the hot component. 
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Figure 8. Best fitting metallicities for the spectral fits from multi- 
temperature (5 X VAPEC) and cooling flow model fits (VAPEC+VMCFLOW 
with free minimum temperature) to the data extracted from the 99 per cent 
PSF. Note that N, O and Fe are allowed to vary between the two components 
in the VAPEC+VMCFLOW model. 
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Figure 9. Best fitting spatial size of spectral fits using the line widths for 
the 5x VAPEC model (adding the PSF contribution). The lines connect tem- 
perature components which have their sizes tied together 
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Figure 10. Comparison of the RGS with Chandra emission measures per 
unit temperature as a function of temperature. The Chandra results were 
extracted from a 2 arcmin radius around the core of the cluster The solid 
line shows the expected distribution of emission measure for a lOM, ;; yr^' 
isobaric cooling flow at Solar metallicity, cooling from a temperature of 
4.5 keV to 0.0808 keV using the MKCFLOW model. 



4.3 Multi-temperature model 

We can extend the two-temperature model with more temperature 
components. Multitemperature fits with several free temperatures 
typically become unstable when spectral fitting. To avoid this we 
constructed a model containing components with fixed tempera- 
tures, allowing the emission measure of each component to vary. 
We used a multi-temperature model containing five components 
(5 X VAPEC). We used a range of temperatures in the model to ac- 
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Figure 11. Maps of the Chandra emission measure per unit area using multi 
temperature model. The circles show the best fitting HEW of the emission 
from the line widths in the spectral fits to the RGS data (the 0.4 and 0.8 keV 
components are forced to be equal, as are the 2.4 and 3.2 keV components), 
including the effect of the PSF. The RGS cross-dispersion direction for the 
longer observation lies along the cool plume. 



count for the ranges of temperatures in the cluster the spectrum is 
sensitive to, with components at 0.4, 0.8, 1.6, 2.4 and 3.2 keV. We 
also tried an 8 component model, but the quality of the fit was only 
slightly improved and it was difficult to constrain the parameters. 

Table|2]show the best fitting parameters for these two models. 
To reduce the number of free parameters, we only used a single set 
of metallicities for the temperature components (tieing O, Ne, Mg, 
Si, Ca, Fe and Ni). We however split the components by tempera- 
ture to allow for two N values, otherwise there are obvious resid- 
uals around the N VII line. Rather than allow the spatial scale of 
each of the components to vary, or fixing them to be all the same, 
we apply three different spatial smoothing scales to groups of the 
temperature components. 

We plot the abundance of each variable relative to solar in 
Fig. [8] for the 5 x VAPEC model. In Fig.|9]is shown the best fitting 
HEW sizes from the line widths against the temperature of the com- 
ponent. We note that velocity broadening of the lines could mimic 
spatial broadening, however, velocities above 700 km s^' are re- 
quired to make a substantial difference to these measurements. 

The emission measure of each component per unit tempera- 
ture (assuming a temperature bin size which is half the logarithmic 
change to the next adjacent temperature value) is plotted as a func- 
tion of temperature in Fig.[T0] In this plot we also show values from 
a Chandra spatially-resolved analysis of the inner 2 arcmin. These 
were created from spectra extracted from contour-binned ( |Sanders| 
2006 1 regions containing around ~ 10^ countQ The spectra were 



' These are the same regions and spectra used to create figures 6 and 7 in 
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fit with a multi-temperature 5xAPEC model with the same five 
fixed temperatures as the RGS fit, but assuming Solar abundance 
ratios. The model assumes each temperature component in individ- 
ual regions has the same metallicity. 

The RGS and Chandra results are roughly similar, considering 
we are comparing the Chandra results in the inner 2 arcmin radius 
with the larger field of view of the RGS instruments and the decline 
of the RGS effective area off-axis. The one major discrepancy is 
the apparent lack of any gas above 3 keV in the RGS results (see 
Section [5!4| for a discussion). We also plot the HEW size from the 
line widths on the Chandra emission measure maps on Fig. [TT] 
showing they are comparable. 

The solid line in Fig. [TO] shows the expected distribution of 
emission measure for a simple isobaric cooling flow, cooling at the 
rate of lOMy, yr^' without any heating. It can be seen that the ob- 
served emission measure of gas decreases more steeply with tem- 
perature than predicted by a simple cooling flow. 



4.4 Single cooling flow component 

We have investigated how well the observed spectrum can be mod- 
elled by a simple cooling flow. Our first model was based on a 
VAPEC component plus a VMCFLOW cooling flow component. We 
used the new feature in XSPEC version 12 to base the cooling flow 
model spectrum on a APEC thermal model rather than a MEKAL 
one. We note, however, that this form of the model is not internally 
self consistent, as the quantity of gas at each temperature is com- 
puted by assuming the luminosities of the MEKAL model. We made 
our own consistent version of the model, but this had no effect on 
the predicted spectrum, so we show results from the XSPEC model 
here. 

We tried two forms of the model: a full cooling flow where 
the lower temperature of the cooling flow was constrained to be 
the minimum possible (0.0808 keV) and the second reduced model 
where it was allowed to be a free parameter. In this model we al- 
lowed the N, O and Fe metallicites to vary between the VAPEC and 
VMCFLOW components, but fixed the other metallicities to have the 
same values in the two components. 

The best fitting parameters for the two models are shown in 
Table |2] The full cooling flow model gives a mass deposition rate 
of 7.4 ± 0.5 Mp, yr^' , whereas the reduced model obtained a rate 
of 8.7 ±0.5 M0 yr"' coohng to 0.54± O.OI keV. The metallicities 
of the cooling flow component were lower for the model where 
the gas cools to the minimum value, presumably to decrease the 
strength of the emission lines. 

The reduced cooling flow model gives a substantially better 
quality of fit to the spectrum than the full model (x^ = 6381 ver- 
sus 6494), but is poorer than the five component multitemperature 
model ix^ = 6362). 
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Figure 12. Mass deposition rates in the absence of heating using a cooling 
flow model made up of five different temperature ranges (5 x VMCFLOW). 



steps in temperature from 3.2 to 2.4, 2.4 to 1.6, 1.6 to 0.8, 0.8 to 
0.4 and 0.4 to 0.0808 keV. The metallicities are assumed to be the 
same in each component, except for the N metallicity in the three 
coolest components, which is allowed to vary separately from the 
hotter components. 

We show in Fig.[T2]the mass deposition rate in the absence of 
heating for each of the cooling flow components. We measure mass 
deposition rates down to 0.4 keV and an upper limit of 0.8 M r, yr^ ' 
below that temperature. The results are consistent with the multi- 
temperature model results. There is systematically less gas detected 
at lower temperature than expected from a simple cooling flow (as 
inFig.flOl. 



4.6 Direct spectral line measurements 

Instead of direct spectral fitting to the whole of the spectrum, the 
strength of the emission lines can be used to gauge the amount of 
cooling taking place through the temperature range they are sensi- 
tive to. The amount of flux in a line can be compared to that ex- 
pected from a cooling flow model. As there are several possible 
lines which can be examined, this allows independent determina- 
tions of the mass deposition rate, but does not use the full spectral 
information available when spectral fitting. 



4.5 Multiple cooling flow components 

A simple cooling flow is unlikely to be a good model to the complex 
distribution of gas in the core of the Centaurus cluster, even in the 
presence of cooling gas. We have therefore used a model where the 
cooling flow is split into different temperature ranges and we allow 
the mass deposition rate to vary in each range. We use fixed ranges 
in temperature to obtain a stable fit. The model examined here has 

[Sanders & Fabian] j2006a) , extracted from Chandra observation IDs 504, 
504,4954,4955 and 5310. 



4.6.1 Fe XVII 

To examine the Fe XVII line powers we fitted a spectral model made 
up of an absorbed powerlaw plus Gaussians, fitted to the first order 
spectra between 14 and 18.5A extracted from 99 per cent of the 
PSF. We used a Gaussian for each of the distinct Fe XVII lines (Ta- 
blellJ, two Gaussians at 16.07 and I6.OOA for Fe XVIII and O VIII, 
and Fe XVIII Gaussians at 14.2lA and 17.62A. The strong lines 
at 15.01, 16.00, 16.07, and 17.06 were allowed to have variable 
widths. The other lines were fixed at zero width. Galactic absorp- 
tion was fixed at 8.56 x 10^" cm^^ . 
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Table 3. Line luminosities for some of the more interesting lines. The listed wavelengths are rest wavelengths, averaging the position of lines in blends. The 
temperature range is the approximate range in temperature of gas the lines are emitted from (where the emissivity is within an order of magnitude of the peak). 
The widths are measurements of the Gaussian width, upper limits (2a) or fixed at 0. The power and mass deposition rates are corrected for Galactic absorption. 
The mass deposition rates are calculated assuming a cooling flow cooling to 0.0808 from 4.5 keV temperature with Solar metallicities. 

Line Wavelength (A) Temperature range (keV) Width (A) Power (10^' erg s^' ) M(M0yr^') 



FeXVII 15.01 0.2^0.9 < 0.07 12.4±1.5 1.6 

Fexvil 15.26 0.2^0.9 3.6±1.3 1.6 

FeXVII 16.78 0.2^0.9 4.3±1.2 1.7 

Fexvil 17.06 0.2^0.9 0.047 ±0.009 19.3 ±2.0 3.0 

Nvil 24.78 <0.85 0.050±0.014 11.5 ±1.6 9.7 

Ovil 21.60 <0.4 <4.6 < 3.4 

Ovil 22.10 <0.4 <7.7 < 2.9 




14.5 15 15.5 16 16.5 17 17.5 18 
Wavelength (A) 

Figure 13. Absorbed powerlaw plus Gaussians model fit to the Fe XVII 
region of the spectrum. The top panel shows the background-subtracted data 
and best fitting model. The bottom panel shows the contribution, Aj^, of 
each channel to x^- The RGSl and RGS2 data and models are shown in 
black and red, respectively. The spectrum has been rebinned to a signal-to- 
noise ratio of 10 in XSPEC. 
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Figure 14. Predicted ratio of the 15.01 line to the blended 17.06 line as a 
function of temperature for a variety of different plasma codes. The line 
with suiTounding shaded region shows the best fitting value and its uncer- 
tainty. 



Table [3] shows the best fitting widths and de-absorbed line 
luminosities. Fig. |13| shows the data with best fitting model. The 
line powers would be increased if there were any internal absorp- 
tion within the cluster (as suggested by Crawford et al.|2 005 and 
[Sanders & Fabian|2006b| l. We convert these powers into mass de- 
position rate by comparing them against the expected flux from 
a cooling flow model. This was done by firstly simulating a high 
signal-to-noise spectrum using a lOM^yr^' VMCFLOW model, 
cooling between 4.5 and 0.0808 keV, with Solar abundance, ab- 
sorbed with Galactic absorption. We then measured the fluxes in 
the simulated lines by fitting the same absorbed powerlaw plus 
Gaussian model as we did to the data. The ratio of the line flux 
in the data compared to the simultated model was used to obtain a 
mass deposition rate. Note that the assumed relative abundance can 
make a substantial difference to these values. The mass deposition 
rates obtained for three of the lines agree at --^ 1.6 Mq yr^ ' , but the 
17.06A blended line is substantially stronger than the others. 

The APEC model predicts that the 15.0lA line should be more 
powerful than the 17. 06 A blended line over a wide range of tem- 
peratures. However in our data the 17.06A lines are 56 per cent 
brighter than the 15.0lA line. A possible cause for this is resonant 



scattering, as the IS.OlA transition is resonant with a large oscil- 
lator strength. If this were the case the scattered radiation would 
have to be absorbed, as scattering only increases the spectral width 
of the line, redistributing the flux. 

However, a likely reason for the discrepancy in brightness may 
be due to uncertainties in spectral models. We plot in Fig. |14|the 
predicted ratio of the 15.0 to blended 17. OA lines versus the tem- 
perature for some different plasma codes, including APEC, MAPEC 
(modified APEC results using a many-body perturbation theory 
method to calculate new values for the emissivities of Fe and Ni 
L-shell lines; 'Gu'2007|l, SPEX (the latest MEKAL results from the 
SPEX package; ,Kaastra & Mewe^2000 ) and CHIANTI ( |Dere et alT] 
[l997||Landietal.|2006| l. ' ' 

It appears APEC cannot reproduce the observed line ratio (as 
found by |Xu et al.|2002| >. The other plasma codes are able to and 
agree fairly well. The line ratio appears to indicate that the Fe XVII 
line emission comes from an average temperature of 0.25i'^Q Qy keV. 
Note that although APEC appears to fit the data worse here, its qual- 
ity of fit to the total spectra is at least ^'x} ^ 300 better than the 
SPEX MEKAL model overall. 
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Figure 15. Flux ratio of 21.60A O VII ratio to 15.0lA Fe XVII line. The 
model lines show the ratio as function of temperature for Solar abundance. 
The horizontal line shows the observed upper limit. 



The Fe XVII I7.06A line width of 0.047A translates into a 
HEW of 36 arcsec (including the effect of the mirrors). This is 
larger than the width of the emission line in the cross-dispersion 
direction (Fig. [5] compatible with the PSF). The longer observa- 
tion analysed here has the cross-dispersion direction placed along 
the plume. The difference in source extent between the line widths 
and cross-dispersion profiles may be due to the very coolest gas be- 
ing extended in the direction perpendicular to the cross-dispersion 
direction (Fig.[TTJ C rawford et al.|2005> . 

If a cooling flow model plus powerlaw continuum is fit to this 
region of the spectra containing the Fe XVII lines, the best fitting 
mass deposition rate is 1.90 ± 0.13 yr^' . This is for a cool- 
ing flow at Solar abundance cooling from 4.5 to 0.0808 keV. As 
the Fe metallicity decreases, the mass deposition rate climbs to 
3.8M0yr-i. 



4.6.2 O VII 

We next placed upper limits on the O VII lines at 21.6 and 22.1 A. 
We fitted an absorbed powerlaw to the spectra from 19.5 to 23A 
plus two zero-width Gaussians at the line positions. We calculated 
upper limits for the line powers and converted these to upper limits 
on the mass deposition rate. These upper limits of ~ SMeyr^' 
are compatible with the Fe XVII results, although Solar metallicity 
is assumed. The 0/Fe metallicity will affect these values as Fe is 
the main coolant in this temperature range ( [Bohringer & Hensler| 
[T989l l. 

We note that if the emission only came from gas at 0.25 keV, 
we would expect the O VII lines to be much stronger than Fe XVII 
(see Fig. [2] lower panel). The relative strength of the Fe XVII and 
O VII will depend on the metallicity and ionisation balance, how- 
i the expected ratio of the 21.60A O VII to 



15 



ever. We show in Fig. 
I5.OIA Fe XVII line as a function of temperature assuming Solar 
abundance. The ratio of any of the O VII to Fe XVII lines is a strong 
function of temperature. For sensible ranges in abundance (O/Fe 
from 0.3 to 1 Solar), it appears unlikely that the temperature of the 
material emitting the Fe XVII lines is less than 10^'^^ - 10^'^ K. 
O/Fe of around O.I Solar is required for 0.25 keV material to be 



Figure 16. Flux ratio of 24.78A N VII ratio to 15.01A Fe XVII line. The 
model lines show the ratio as function of temperature for Solar abundance. 
The shaded region shows the observed ratio. The vertical line shows the 
minimum temperature from the O VII to Fe XVII Une ratios (Fig.jTs}. 



consistent with the line ratio. This means a likely lower temperature 
is 0.3 — 0.35 keV, still within the allowable range of the Fe XVII ra- 
tio. These estimates all assume that the line emission comes from 
isothermal gas. 

We also note that resonant scattering could be important for 
the 21.6A Ovil line, which has a high oscillator strength (0.7). Al- 
though scattering by itself would not decrease the line flux, it would 
broaden the spectral line, making it harder to detect. Absorbing ma- 
terial within the scattering region could also reduce the line flux. 



4.6.3 N VII 

We also try to measure the mass deposition rate in the absence of 
heating using the N VII line. We fit the region 23 to 30A with an 
absorbed powerlaw plus Gaussian with variable width. The power 
in this line compared to the Fe XVII lines is large: assuming so- 
lar abundance its power corresponds to a 10 M0 yr^' cooling flow. 
This indicates that the N abundance of the cool X-ray emitting gas 
relative to Fe is large, as was also shown by the spectral fits. Addi- 
tionally the N VII line covers a larger range in temperature, and so 
is sensitive to hotter gas which shows a higher effective mass de- 
position rate (Fig.|12[). We note that any additional absorption will 
increase this line power relative to the Fe XVII lines. 

The spectral width of the line is also very close to the Fe XVII 
17.06A line width, meaning it originates from the same spatial re- 
gion. The finite widths of the lines (a point source would provide an 
upper limit as the response already takes account of the XMM PSF) 
show that the emitting region is resolved, with a HEW of around 
35 arcsec (taking into account the XMM PSF HEW). 

If the lower bound on the temperature is 10^'^^ K, then the 
ratio of the N XVII to Fe XVII 15.0lA line implies supersolar N 
abundances (Fig. |16^ . The minimum N/Fe metallicity as indicated 
by the APEC and SPEX model line ratios is at least 4Z0. This value 
is in rough agreement with the spectral fitting results. 
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5 DISCUSSION 

5.1 Temperature range 

The observed range in temperature of X-ray emitting gas, from 
3.7 keV ( [Sanders & Fabian|2002j to ~ 0.35 keV is one of the largest 
observed up to now in a cluster of galaxies. These RGS results con- 
firm earlier CCD detections of cool gas. 

Another cluster with a large apparent X-ray emitting temper- 
ature range is the Perseus cluster, containing X-ray emitting fila- 
ments at around 0.7 keV in a cluster up to 7 keV ( [Sanders & Fabian] 
[2007 | l, although these are clearly influenced by magnetic fields. 

The result emphasises that gas certainly does exist at temper- 
atures below the often quoted factor of 2 or 3 in galaxy clusters. 
Fig. [To] shows that the emission measure of cooler components is 
progressively smaller than a cooling flow model as the temperature 
decreases. The lack of detections of cool gas in other clusters using 
RGS observations could be due to their relatively short exposures, 
which are typically only 30-40 ks ( [Peterson et al.|2003) . 

The large range in temperature in this cluster may be con- 
nected to the lack of disruption in its core. The high enrichment 
of the central parts of the cluster indicates it has been stable for 
8 Gyr or more ( [Sanders & Fabian[2006a) . 

5.2 Cooling flow model 

Fig.[TO]and Fig.[T2]show that a full simple cooling flow, where gas 
is cooling by radiative cooling at a constant rate, is not operating 
in this cluster. However there are some possible reasons why the 
amount of cold gas may be underestimated: 

(i) These models assume that there is a fixed metallicity of mate- 
rial as a function of temperature. The O and Fe abundances change 
the cooling function and strengths of lines significantly. Fig. [TT] 
shows how the mass deposition rate for the 5 x VMCFLOW cooling 
flow model varies as a function of metallicity, for the two coolest 
temperature ranges. This was made by iterating over a grid of O 
and Fe abundances and measuring the mass deposition rate at each 
point. It shows that the real cooling rate is a strong function of 
metallicity. If the metallicity of the material in the central regions is 
small there could be significantly more cooling taking place. There 
is evidence for a drop in metallicity in the central regions jSanders[ 
[& Fabian[2002[ l. This analysis highlights that the best fitting Fe and 
O metallicities for the two coolest components are small. 

Indeed, O VII emission is the main indicator of cool gas in 
this waveband, but Fe line emission is actually the main coolant 
( [Bohringer & Hensler|1989} . Therefore the O lines are suppressed 
if the Fe metallicity is increased while the O metallicity is constant, 
for a fixed observed Fe XVII strength in a cooling flow. Therefore 
gas below 0.35 keV can exist without strong O VII emission. 

(ii) Intrinsic absorption could hide much of the cooling taking 
place. There is evidence that there are large amounts of obscuring 
material in the core of this cluster ( [Crawford et al.|2005| [Sanders[ 
[& Fabian|2006b l, where the coolest material resides. It is difficult 
for spectral fitting to account for this as the continuum is hard to 
detect for low-temperature gas. We investigated this by fitting the 
5 X VMCFLOW to the data with varying levels of absorption applied 
to the two lowest temperature components (note that these results 
were obtained by varying the absorption separately for each com- 
ponent). In Fig.[T8]we plot the increase in mass deposition rate it is 
possible to obtain by increasing the absorption. 

(iii) The cooling flow model fitted to the data assumes that iso- 
baric (constant pressure) cooling is taking place. If the cooling time 




Iron metallicity (Solar) 



Figure 17. Variation of mass deposition rate with Iron and Oxygen metal- 
licity for the two coolest temperature ranges in Fig. |12| The top panel shows 
as contour lines the best titling mass deposition rates in M, : yr"' for the 0.8 
to 0.4 keV temperature range. The bottom panel shows upper limits for 0.4 
to 0.08 keV. Also plotted as shaded regions on both panels are the 1, 2 and 
3-fT confidence regions for the O and Fe metallicity, calculated using the 
method from the minimum value. 
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Figure 18. Variation of mass deposition rate with excess column density 
above 8 x 10-" cm^^. Shown are the best fitting mass deposition rates in 
M .) yr^' for the 0.8 to 0.4 keV temperature range, or upper limits for the 
0.4 to 0.08 keV temperature range. 
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becomes shorter than the sound-crossing time, then isochoric (con- 
stant density) coohng may be more appropriate. This could also 
result from the magnetic pressure in clumps becoming dominant 
as the cooling proceeds. The strength of UV emission lines from 
isochoric cooling is only around 60 per cent of that expected from 
isobaric cooling ( [Edgar & Chevalier|1986| l. This would lead to the 
amount of cooling being underestimated if a cooling flow were tak- 
ing place. As we show below, the mean radiative cooling time is 
only ~ lO' yr for the very coolest X-ray detected gas. 

(iv) If the X-ray emitting gas mixes with the cooler optical line 
emitting material in the very central regions ^Crawfordetal.|2005| l, 
than it will rapidly cool non-radiatively much faster than a cool- 
ing flow model would predict ( [Fabian e t al. 2002). The thermal 
energy is eventually radiated at much longer wavelengths, in the 
optical/UV/IR filaments, for example, or by dust. There is, indeed, 
extended mid-IR emission observed from NGC 4696 jKaneda etal.[ 
[2007 [ l, which could in part be powered by such a mixing process. 

When fitting the data with a cooling flow model, we find pro- 
gressively smaller mass deposition rates at lower temperatures, as 
found previously by Peterson et al. ( 2003 ). However we detect gas 
over a much wider range in temperature (a factor of 10 rather than 
3) than in their shorter observations, with a wider range of best fit- 
ting mass deposition rates (a factor of 50 rather than 10), in reality 
corresponding to detecting a wider range of emission measures for 
the temperature components. 

Taking gas at 0.35 keV in pressu re equilibrium with the ga s 
with «e ~ 0.1 cm"^ and kT ~ 1 keV JSanders & Fabian|2006b|, 
implies a density of this cool material of around «e ~ 0.3 cm^ . 
We can calculate the bolometric luminosity from a unit volume of 
this material (using the abundances from the 5 x VAPEC model) and 
its internal energy (assuming |fcB^"tot)- Dividing the two gives 
a mean radiative cooling time of only ~ 10^ yr for this 0.35 keV 
material. 

The 0.4 keV component in the 5 x VAPEC model is close to 
0.35 keV in temperature. Assuming the material is in pressure equi- 
librium with its surroundings and taking its emission measure from 
the 5 X VAPEC spectral fit, we can estimate it occupies a volume 
of 2.3 kpc-'. The region that the Fe XVII lines are emitted from 
has a radius of at most 10 arcsec, implying a maximum volume 
of 120 kpc^ This means the 0.4 keV material must fill only 2 per- 
cent of that very central region. Chandra temperature maps of the 
core show that the coolest gas is not volume-filling, but exists in the 
form of cool blobs ( [Crawford et al.|2005| l. 



5.3 Heating 

The amount of heating which would be required to offset cooling of 
the ICM in the cluster can be estimated. We took the 5 x VMCFLOW 
cooling flow model fits to the spectra and generated a chain of pa- 
rameters using a Markov Chain Monte Carlo. The chain was cre- 
ated with the built-in XSPEC functionality and was checked for con- 
vergence by repeating the analysis with a different starting point. 
We generated luminosities for each of the sets of parameters in the 
chain for the individual temperature ranges in the model. We ex- 
amined the distribution of luminosities to estimate the medians and 
uncertainties. The luminosities of each temperature range are the 
rates of heating required to offset cooling. 

In Fig. [19] we show the amount of heating required for each 
temperature range of gas in order to combat cooling. At tempera- 
tures of above ~ 1.5 keV, significant heating rates of a few times 
lO'*^ erg s^^ are required. For the coolest gas it is likely that mix- 
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Figure 19. Heating rate required to combat cooling, calculated for each of 
the individual temperature ranges in the 5 x VMCFLOW model (see Fig. [12^ . 

ing and conduction complicate this picture. The cool gas we ob- 
serve here must exist as distributed clumps and the short radiative 
cooling time places constraints on heating models. 

5.4 Metallicity of the gas 

The best fitting metallicities we measure show some model depen- 
dence. Those models with larger number of temperature compo- 
nents tend to have higher metallicites. However there is quite good 
agreement between the ratios of the various elements between the 
best fitting models here and with the previous CCD results. 

The metallicities we measure here are smaller than the peak 
values obtained using CCD data from Chandra and XMM ( [Sanders[ 
[& Fabian[20 06a'). It is however difficult to compare them exactly as 
the CCD results are spatially-resolved and the metallicity appears 
to rise to high values, but drops in the very central region. The 
temperature is also declining towards the centre. Over the region 
examined by the RGS detectors, it is likely that there is a rise in 
metallicity with radius. 

Ni shows significantly smaller metallicities of around 2.4 Zq 
here when compared to earlier Chandra and XMM spatially- 
resolved CCD results (which peak near 4 Zq; 'Sand ers & Fabian[ 
|2006al. This can be mostly resolved by allowing the Ni to vary 
between the coolest and hotter components, obtaining values of 
3.5+Q4 and I.T+q jZo, respectively (for the 5xVAPEC model). 
This highlights that the assumptions or choices we made in mod- 
elling the data can have a large impact on the obtained quantities. 

The generally smaller metallicites we find here, when com- 
pared to CCD results, may be due to problems modelling the con- 
tinuum from hot gas outside of the central region. Fig. [lO] shows 
that the RGS results require very little hot gas in the spectral fits 
when compared to Chandra. If we add the expected 3-4 keV hot 
component from the Chandra spectral fits to the RGS spectral mod- 
els, freezing its emission measure to the same value, then we typ- 
ically find the metallicities we measure are increased by 30 to 
40 percent. The RGS spectra appear to require much less of this 
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Figure 20. Ratio of metallicity of elements to Fe obtained from the RGS 
data witli tlie SxVAPEC model, compared to the results from the central 
regions of the cluster from combined Chandra and XMM EPIC MOS data 
tSanders & Fabian|2006a) . 



hot material than all previous analyses. Therefore there may be cal- 
ibration uncertainties at short wavelengths in the RGS spectra. 

Another possible reason for higher continuum, leading to un- 
derestimates of the metallicity, would be the previously observed 
non-thermal emission ( [Alien et al. 2000, Di Matteo et al. 20001. 
The sedimentation of helium in the core of the cluster would also 
lead to higher continuum l |Ettori & Fabian|2006| l. 

The ratio of elemental abundances is likely to be more robust 
than measurements of individual abundances. In Fig. |20|we plot 
the ratio of the metallicity of various elements to Fe (the ratios are 
Zq/Zq), measured using the multitemperature SxVAPEC model 
(Fig. [Sj. Also shown are the same ratios obtained from Chandra 
and XMM EPIC MOS data. These were calculated by averaging 
the central XMM and Chandra East and West points from figure 15 
(bottom panel) in [ Sanders & Fabiari) (2006a ). The plot shows that 
we get good agreement between the RGS, and previous CCD ratios 
within the inner 10 kpc radius, with some discrepancy for Si/Fe 
and Ni/Fe. These results indicate that the CCD-derived metallicity 
ratios in' Sanders & FabiM i'(2006al are robust. 

If the coolest gas is examined separately from the hotter gas, 
it appears to have lower Fe and O metallicity, confirming the cen- 
tral abundance drop in the spatially resolved CCD measurements. 
Fig. [it] shows the best-fitting values of the cooler gas are only 
around 60 per cent of the metallicity of the hotter gas. 

It is possible that the metallicity variations are with tempera- 
ture, not radius, although we do not understand why. If there is a 
continued decrease of metallicity as temperature decreases then the 
low O VII result might be understood. It is unlikely that the low 
metallicity gas can be introduced into the centre, and we presume 
that the metals are removed. At a temperature of 0.2 keV the cool- 
ing rate is sensitive to the metallicity. One simple explanation is 
that the high metallicity gas is inhomogeneous, conduction is sup- 
pressed between clumps, leading to the the high metallicity regions 
cooling out ( [Morris & Fa bian 20031. This means that both the in- 
jection of metals is a inhomogeneous process and the metals remain 
poorly mixed for a considerable time (several Gyr). 

We measure values for the N abundance of between 1.6 (from 
the 5 X VAPEC model) to 4 Zq (for the 5 x VMCFLOW and lower limit 
from the line ratios). 

The nitrogen emission lines in the optical spectrum of the 
emission-line filamentary nebula in NGC 4696, the central galaxy 



in the Centaurus cluster also indicate a high nitrogen abundance. 
Specifically, the [NII]/Ho! line ratio found by [Johnstone et al.[ 
( [1987^ of 3.2 is higher than the ratio of less than two which can 
be accounted for by photoionization or shock models for a solar 
abundance plasma (Heckman et al.|[1989[ [Voit & Donahue|1990[ 
[Crawford & Fabian|l992^ high nitrogen abundance in both the 
optical filaments and the hot gas supports the hypothesis that they 
have a common origin, with the filaments representing gas which 
has cooled from the hotter phase. 



6 CONCLUSIONS 

We clearly detect Fe XVII emission from the core of the Centaurus 
cluster Fitting spectral models and measuring emission line ratios 
shows that the temperature of the gas declines to 0.3 to 0.45 keV. 
These results from this deep RGS observation show the widest 
range in ICM temperature detected in a cluster of galaxies, with 
a factor of 10 or more. 

The data confirm that the metallicity of the ICM declines in 
the very central regions. Whether this is a real decline in metallicity, 
the result of an inhomogeneous metallicity distribution ^Morris &] 
[Fabian[2003 l) or excess continuum is unclear We do, however, find 
strong nitrogen enhancement in the inner 6 kpc radius. 

The very coolest gas is concentrated in the centre of the clus- 
ter, as previously found by Chandra, but there is much less cool 
material than would be expected from a simple cooling flow. 

The problem of how to prevent this very cool gas, which has a 
mean radiative cooling time of 10^ yr, from cooling remains. Any 
mechanism which heats the gas must preserve the sharp metal- 
licity gradients and not switch off for timescales of greater than 
10^ yr Alternatively, the ICM in the cluster may be cooling in a 
non-radiative way below 0.4 keV. If this material formed stars, the 
rate of mass deposition is compatible with the metallicity of the 
ICM (San ders & Fabian|2006at . 

Using this very deep observation we are able to measure the 
amount of cool X-ray emitting gas very precisely in a cluster core. 
In order to get a more complete picture, it is necessary to make 
deep observations of several nearby clusters to see whether the 
large ranges in temperature are found in all objects. 
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